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ABSTRACT 

Starting from the ^50000 quasars of the Sloan Digital Sky Survey for which Mgn line 
width and 3000 A monochromatic flux are available, we aim to study the dependence 
of the mass of active black holes on redshift. We focus on the observed distribution 
in the FWHM-nuclear luminosity plane, which can be reproduced at all redshifts 
assuming a limiting Mbh, a maximum Eddington ratio and a minimum luminosity 
(due to the survey flux limit). We study the z-dependence of the best fit parameters 
of assumed distributions at increasing redshift and find that the maximum mass of 
the quasar population evolves as log(M BH ( max )/M0) ~ 0.3z + 9, while the maximum 
Eddington ratio (<~ 0.45) is practically independent of cosmic time. These results are 
unaffected by the Malmquist bias. 
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1 INTRODUCTION 

, In the last years a substantial effort has been devoted to 
measure black hole (BH) masses for various quasar samples 
covering a wide range of redshifts and luminosities. McLure 
& Dunlop (2004), from H/3 and Mgn, measured virial black 
hole masses (Mbh) for ~ 10000 quasars with z ^ 2.1 in- 
cluded in the Sloan Digital Sky Survey (SDSS) Data Re- 
lease 1 (DR1). Fine et al. (2006) used composite spectra to 
measure the dependence on redshift of the mean BH mass 

, for an L* subsample of the 2QZ quasar catalogue (Croom 
et al. 2004) from z ~ 0.5 to z ~ 2.5. Shen et al. (2008) 
listed BH masses for ~ 60000 quasars in the redshift range 
0.1 < z < 4-5 contained in the SDSS DR5, by means of virial 
BH mass estimators based on the H/3, Mgn and Civ lines. 

A common result of these works is that the mean BH 
mass of the QSO population at given z appears to increase 
with redshift, but the observed z-dependence is dominated 
by the well-known Malmquist bias, because the BH mass 
strongly correlates with the central source luminosity (see 
Vestergaard et al. 2008 for a detailed analysis of the selection 
bias effects). McLure & Dunlop (2004), for instance, suggest 
that the observed active BH mass evolution is entirely due 
to the effective flux limit of the sample. 

A full understanding of this scenario would give impor- 
tant insights on the BH formation and evolution and on the 
activation of the quasar phenomenon. Moreover, along with 
a parallel study on the dependence on redshift of the host 
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galaxy luminosity (mass), this would enlighten on the joint 
evolution of galaxy bulges and their central black holes. For 
these reasons, it is of focal importance to trace the depen- 
dence of Mbh on z, overcoming the problems related to the 
Malmquist bias. 

We start from the recently published SDSS DR5 quasar 
catalog (Schneider et al. 2007) and focus on the ~ 50000 
quasars for which Mgn line width and 3000 A flux are avail- 
able (Shen et al. 2008). The sample selection is described 
in Section [5] The sample (0.35 < z < 2.25) is divided in 8 
redshift bins and it is shown that, in each bin, the object 
distribution in the FWHM-luminosity plane can be repro- 
duced assuming a minimum luminosity, a maximum mass 
and a maximum Eddington ratio (Sections 13.21 and 13. 3[) . 
Comparing the assumed probability density to the observed 
distribution of objects, the parameters can be determined 
in each redshift bin (Section 13. 4[) . This procedure is shown 
to be unaffected by the Malmquist bias (Section 13. 5p . and 
provides a method to study the "unbiassed" dependence on 
redshift of quasar BH masses and Eddington ratios (Section 
14- 1 p - In Section ^. 2l we test the dependence of our results on 
the assumed tblr — ALa calibration. We compare our results 
with previous literature in Section \4. 31 and in Section ^, 4l we 
discuss some implications of our findings for the study of the 
co-evolution of supermassive BHs and their host galaxies. A 
summary of the paper is given in the last Section. 

Throughout this paper, we adopt a concordant cosmol- 
ogy with H = 70 km s" 1 Mpc -1 , fi m = 0.3 and fi A = 0.7. 
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Figure 1. Average values of the absolute magnitude Mj (trian- 
gles) and of Mbh (circles) vs. redshift, in bins Az = 0.15. The 
typical standard deviation of each bin is given in the inset. 



2 THE Mgll SAMPLE 

The SDSS DR5 quasar catalogue (Schneider et al. 2007) 
contains more than 77, 000 quasars. It covers about 8000 
deg 2 and selects objects with Mi < —22, have at least one 
emission line with FWHM larger than 1000 km/s or are un- 
ambiguously broad absorption line objects, are fainter than 
i = 15.0 and have highly reliable redshifts. 

Shen et al. (2008) calculated BH masses for ~ 60000 
quasars in the redshift range 0.1 < z < 4.5 included in 
the SDSS DR5 quasar catalogue, using virial BH mass esti- 
mators based on the H/3, Mgn and Civ emission lines. They 
provide rest-frame line widths and monochromatic luminosi- 
ties at 5100 A, 3000 A and 1350 A (see Shen et al. 2008 for 
details on calibrations, measure procedures, corrections). 

In the following we will focus on the ~ 50000 quasars 
from the Shen et al. (2008) sample for which Mgll line width 
and 3000 A monochromatic flux are available (hereafter Mgll 
sample). We assume the virial theorem and adopt the cal- 
ibration of McLure & Dunlop (2004) to evaluate the BH 



log M BH = 6 + log(a) + 21og(FWHM) + b log ALa 



(1) 



with a = 3.2±1.1 and b = 0.62±0.14. Here Mbh is expressed 
in solar masses, FWHM in units of 1000 km/s and ALa in 
units of 10 44 erg/s. 



3 DESCRIPTION OF THE PROCEDURE 
3.1 Malmquist bias 

Figure Q] shows the mean absolute i magnitude (see Shen 
et al. 2008 for details) vs. redshift of the Mgn sample. The 
effects of the Malmquist bias are apparent as an increase 
of the average observed luminosity with redshift. The mean 
BH mass vs. redshift is overplotted to the mean Mi(z): it 
is apparent that the average observed BH masses follow the 
same trend as the absolute magnitudes with a higher dis- 
persion, as expected given that the distribution of the line 
widths does not depend on luminosity or redshift (Shen et al. 



2008). This suggests that the ^-dependence of the observed 
BH masses is strongly subject to a Malmquist-type bias, be- 
cause at high redshift one cannot observe low mass objects. 
In order to trace the "unbiassed" dependence of active BH 
masses with redshift, one should consider a combination of 
two effects, namely the z-dependence of the quasar number 
density and the increase of the average mass of quasar popu- 
lations with redshift. To illustrate these effects consider two 
extreme cases: 

(i) The Mbh distribution does not depend on redshift, 
but the quasar number density increases until z ~ 2 — 2.5. 
At any redshift there is a population of low mass (~ 10 s Mq) 
quasars, which cannot be observed at high redshift, and the 
population of high mass (~ 1O 9 ' 5 M0) active BHs that is 
observed at z > 1.5 is the high mass end of the Mbh distri- 
bution. 

(ii) The quasar Mbh distribution shifts toward higher 
masses at increasing redshift. The population of low mass 
(~ 10 s Mq) objects that is observed at low redshift is not 
present at all at z > 1.5. The observed increase of Mbh 
with redshift, in active BHs, is "true" and it is not due to a 
Malmquist-type bias. 

Of course, each of these pictures is per se unlikely: the 
observed dependence on redshift of quasar BH masses is 
due to a combination of both these effects. In the follow- 
ing, we will concentrate on these points using statistical 
arguments, focusing on the distribution of objects in the 
FWHM-luminosity plane. 



3.2 Quasar distribution in the FWHM— ALa plane 

Figure [5] shows the objects of the Mgll sample in the 
logFWHM-logAZ/A plane (see Fine et al. 2008 for a simi- 
lar approach). The sample has been divided in 8 redshift 
bins of equal co-moving volume. In each panel, it is appar- 
ent that the data-points form a sort of "triangle", the left 
side of which represents a cut due to the the survey flux 
limit (which gives raise to the Malmquist bias). From Eq.Q] 
the loci of quasars with constant mass are represented in 
this plane by straight lines with fixed slope, as plotted in 
the figure: 



log(FWHM) = -0.31 log ALa + 0.5 log 



M BH 



3.25 



(2) 



where units are the same as in Eq. [T] We propose that the 
top right side of the triangle is representative of a maximum 
mass in the quasar sample. 

The third (i.e. the bottom right) side of the triangle is 
supposedly due to the Eddington limit, as the loci of quasars 
with constant Eddington ratios are again straight lines. The 
dependence of FWHM on the monochromatic luminosity at 
a given Eddington ratio is fixed assuming the bolometric 
correction by Richards et al. (2006; BC3000 = 5.15) and 
Eq. Q] This yields: 



log(FWHM) = 0.19 log ALa 



0.5 log + 0.05 (3) 

LEdd 



where units are the same as in Eq. [T] Note that in each 
redshift bin, the plotted cuts describe qualitatively well the 
shape of the quasar distribution in the FWHM-luminosity 
plane. 
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Figure 2. The 8 panels show the Mgn sample in the FWHM-luminosity plane at increasing redshift. Dotted, dashed and dash-dotted 
lines (and lines parallel to them) represent the loci of constant monochromatic luminosity, constant mass and constant Eddington ratio 
respectively. 
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Figure 3. The shape of P((Z) (Eq.[5J, P m (m) (Eq.[B) and P e (e) 
(Eq.0. 



3.3 Construction of a probability density 

Now we aim to construct a probability density of quasars 
as a function of FWHM and luminosity with a main crite- 
rion of simplicity. We propose that in each redshift bin the 
object density is only constrained by a maximum mass, a 
maximum Eddington ratio and a minimum luminosity due 



to the instrumental flux limit. We then assume a probability 
density of form: 



Pi,FWHM(i, FWHM) = k-Pi(l) ■ P m (m) ■ P e (e) 



(4) 



where k is a normalization constant and each P is assumed 
to be a smoothed step function, which increases from to 1 
(or vice versa) in a range of width a around a fixed value of 
the independent variable. In the following, we describe our 
results assuming P of form (see Figure [3} : 



m) = 



— ( 



exp 

— oo 



{I' — imin — 2(7() 2 



2af 



dl' 



exp 



exp 



Pe(e) = - 

I 

where: 

/ = log XL\ 

m = m(l, FWHM) = lo, 
e = e{l, FWHM) = log ■ 



_ (m — m m ax 
{c Cm ax 



+ 2a m ) 2 



2al 
+ 2a e ) 2 



de 



Lbol 



(•») 

dm' (6) 
(7) 

(8) 
(9) 

(10) 



Here, the parameters Z m i n and a t , m max and a m , e max and a e 
are the minimum luminosity, the maximum mass, the max- 
imum Eddington ratio and the widths of the corresponding 
distributions. These parameters will be determined in the 
following via a best fit procedure. 
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Note that, since the integrals of the P functions diverge, 
we must restrict their domain to use them as probability 
densities (e.g., for values of the parameters I < l m i n + 3ai, 
m > m max — 3<r m and e > e max — 3a e ). This doesn't sig- 
nificantly affect the results, because mass, Eddington ra- 
tio and luminosity are not independent variables (for ex- 
ample, low mass objects also have low luminosities or high 
Eddington ratios), hence the derived probability density 
Pi,fwum{1, FWHM) (Eq. H} is essentially insensitive to the 
shape of P m (m) at low masses, of P e (e) at low Eddington 
ratios or to the shape of Pi (I) at high luminosities. 

3.4 Best fit procedure 

The assumed probability density depends on 6 free param- 
eters, i.e. the minimum luminosity (l m - m ), the maximum 
mass (m max ) and Eddington ratio (e max ) and the widths 
of the corresponding distributions (<?;, a m and <j e ). We fo- 
cus on the first redshift bin and determine the free pa- 
rameters matching with the observed distribution of ob- 
jects in the FWHM-luminosity plane. In detail, for each 
choice in the 6-dimension parameter space, the probabil- 
ity density has been constructed, discretized in boxes with 
A logFWHM=0.04dex and A log AL A =0.2dex and then nor- 
malized to the total number of observed objects, in order 
to evaluate the expected number of objects in each box 
(AlogAZ/A, AlogFWHM). We assumed a Poissonian error 
(i.e. y/n) on the observed number of objects in each box. 
For each choice of the parameters, the expected distribution 
was compared to the observed distribution in the discrete 
log \L\ — logFWHM plane, evaluating the relative \ 2 value. 
The minimum \ 2 determines the best fit parameters. 

In order to determine the uncertainties on these values, 
the same fit procedure was repeated many times compar- 
ing the observed distribution to a set of simulated distri- 
butions of objects, constructed through the Monte Carlo 
method. This procedure allows an estimate of the error 
since we sound the underlying probability density only 
through a finite number of observed objects, the distribu- 
tion of which in the FWHM-luminosity plane ideally fol- 
lows Eq. 4 with a certain random dispersion. In detail, 
given a set of values of the six parameters, we generated 
10 7 points (logALA, logFWHM) with uniform probability 
densities, and then rejected points accordingly to the as- 
sumed ft,FWHM(i,FWHM) at g iven /min, ^^ma X , 6 max , (7;, 
a m and a e , so that the number of simulated points matches 
the number of observed objects. We then calculated the root 
mean square (rms) between the observed and the simulated 
distributions. This operation was repeated for all the possi- 
ble combinations of the six parameters (in a reduced phase 
space around the best fit values). The sextuple which led 
to the minimum rms gave the so-called Monte Carlo best 
fit parameters. This procedure was repeated a dozen times, 
giving as many Monte Carlo best fit values for each param- 
eter, slightly different from one another, but fully consistent 
with the previous determination. For each parameter, the 
standard deviation of this set of best fit values was assumed 
as an estimate of its uncertainty. This uncertainty is much 
larger than that corresponding to A\ 2 = 1. 

In the first panel of Fig. [3] we compare the observed 
distribution with one simulated best fit distribution for the 
lowest redshift bin. It is apparent that the choice of three 



simple distributions in luminosity, mass and Eddington ra- 
tios describes rather closely the data, giving circumstantial 
support to the validity of the virial hypothesis on which the 
theoretical assumptions are based. 

Table [1] (first line) contains the best fit values of the 
6 parameters with relative errors and the reduced \ 2 value 
(hereafter, xt)- The fact that the xt is larger than 1 is inter- 
preted as due to the choice of an oversimplified distribution. 
This doesn't influence our results, because our goal is to 
find a good way to quantify a parameter related to the BH 
mass (and one related to the Eddington ratio) such that it 
is not affected by a Malmquist-type bias (i.e. disentangled 
of the z-dependence of the luminosity instrumental limit). 
In fact, by construction, m max and e max depend neither on 
the quasar number density nor on the survey flux limit (see 
next Section for tests on this statement). 

3.5 Bias analysis and robustness of the procedure 

The effect of the luminosity cut on the results for m max and 
e max can be further tested by simulation. In order to show 
that our results are not affected by the instrumental flux 
limit of the dataset, we selected a subsample from the lowest 
redshift bin applying the probability function Pi (I) (Eq. [SJ 
with a higher luminosity cut, i.e. assuming the Z m i n and o\ 
derived for the third redshift bin in which < z >= 1.22 
(see next Section). This subsample consists of about 1/12 
of the objects in the original lowest redshift sample. The 
fit procedure presented in this paper has been performed 
again on this subsample. Figure [5] (left panel) shows that the 
luminosity cut of a higher redshift bin has negligible effects 
on the results, being these values (m max = 9.20 and e max = 
—0.36) consistent within la with the best fit parameters 
derived for the whole sample (m max = 9.18±0.05 and e max = 
-0.35 ±0.02). 

A similar test has been performed to show that m max 
and e max do not depend on the quasar number density: we 
re-sampled from the first redshift bin rejecting randomly 2/3 
of the objects, in order to obtain a smaller sample with the 
same distribution. The fit procedure was then performed on 
the reduced sample. Again, no significant deviation in the 
determination of the best fit parameters was observed (see 
Figure [5] right panel) . Again, the derived values (m max = 
9.20 and e max = —0.34) are consistent within la with the 
best fit parameters obtained for the whole sample. These 
tests show that m max and e max are independent of the quasar 
number density and of the survey flux limits, and therefore 
indicate that our procedure is not affected by a Malmquist- 
type bias. 

4 EVOLUTION OF THE QSO POPULATION 

4.1 Quasar BH mass and Eddington ratio 
dependence on redshift 

The fit procedure described above is applied to all the red- 
shift bins, in order to determine the best fit parameters and 
their uncertainties as a function of redshift. In each redshift 
bin we compared the best fit minimum luminosity with the 
values inferred through the z-dependence of the luminosity 
distance (see Fig. [6}. It is apparent that the agreement is 
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Figure 4. The 8 panels show the Mgn sample in the FWHM— luminosity plane at increasing redshift: solid black contour plot (levels: 20, 
90, 250 objects per box, see text) represents the discrete observed distribution of objects. Dotted red contour plot (same levels) shows 
the discrete distribution of a sample of objects simulated with the Monte Carlo method, adopting the assumed P; fwhmCi FWHM) 
probability density with the best fit parameters. Dotted, dashed and dash-dotted lines represent Z m i n , m max and e ma x respectively. 



Table X. Best fit values of minimum luminosity, maximum mass, maximum Eddington ratio and widths of the corresponding distributions, 
with errors and xi- The number of degrees of freedom is v = 594 in the first redshift bin (600 data-points and 6 free parameters) and 
v = 595 in the others (600 data-points and 5 free parameters). Data that come from a best fit procedure are displayed in boldface. 



Bin 


< z > 


^min 




'"-max 


o" m 


6max 




xl 


1st 


0.62 


44.30±0.025 


0.26±0.008 


9.18±0.05 


0.31±0.003 


-0.35±0.02 


0.22±0.003 


3.51 


2nd 


0.97 


44.83 


0.23±0.007 


9.35±0.04 


0.31±0.006 


-0.34±0.02 


0.23±0.003 


4.97 


3rd 


1.22 


45.07 


0.23±0.005 


9.42±0.05 


0.31±0.004 


-0.33±0.02 


0.23±0.002 


4.63 


1th 


1.42 


45.24 


0.23±0.008 


9.43±0.05 


0.32±0.003 


-0.34±0.02 


0.22±0.001 


5.11 


5th 


1.61 


45.36 


0.24±0.006 


9.52±0.04 


0.31±0.003 


-0.35±0.01 


0.22±0.003 


4.36 


6th 


1.80 


45.51 


0.22±0.005 


9.60±0.05 


0.32±0.002 


-0.34±0.02 


0.22±0.004 


8.44 


7th 


1.98 


45.61 


0.22±0.005 


9.67±0.05 


0.31±0.003 


-0.33±0.02 


0.22±0.004 


6.84 


8th 


2.15 


45.67 


0.24±0.006 


10.02±0.05 


0.30±0.005 


-0.34±0.02 


0.21±0.003 


7.21 



very good: apart from the highest redshift bin, where the 
3000 A continuum is very close to the red edge of the ob- 
served spectral range and the flux calibration may be unreli- 
able, all the data are consistent with the expectations within 
la. This gives further support to the assumed description 
of the object distributions in the FWHM-luminosity panels 
and suggests to repeat the entire procedure assuming that 
the value of / m in(z) is constrained by cosmology. 

The same fit procedure is then applied again to all the 
redshift bins, but now the dependence on redshift of the 
minimum luminosity is set by cosmology and Z m i n is no more 
treated as a free parameter. In each bin, the Xv was evalu- 



ated and normalized to the xto value obtained in the first 
redshift bin, in order to compare the adequacy of the best 
fit function in the 8 panels. Figure [7] shows that these values 
are almost constant in each redshift bin. Again, Fig. [4] and 
Table [TJ show respectively the Monte Carlo simulated distri- 
butions compared to the observed distributions of quasars 
and the best fit values, their errors and relative \ 2 values. 

The maximum mass and Eddington ratio values are 
plotted versus redshift in Figure [S] Note that the proposed 
Mbh 2-dependence refers to the active BH population and 
not to the total supermassive BH mass distribution. Of 
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Figure 5. The same of Figure [4] referred to a subsample of the 
lowest redshift bin objects, selected accordingly to the luminosity 
cut function of a higher redshift bin (left panel) or randomly 
selected with p = 0.3 (right panel). For comparison, green thin 
lines represent the whole low redshift sample and its best fit Z m i n , 
m m ax and e max lines. 
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Figure 7. Normalized \i values of the best fit function. Filled 
circles refer to the fit procedure described in this work. Open tri- 
angles show the xt values that would be obtained assuming that 
the quasar BH mass is constant with redshift and open squares 
show the xi values assuming that the BH mass evolves with red- 
shift as proposed by Fine et al. (2006). 
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Figure 6. Black dots represent the best fit parameters l m i n vs. 
z, compared to the values expected from cosmology (solid line). 

course, the average mass of the inactive BH population must 
decrease with increasing redshift. 

A linear fit to the maximum mass values (excluded the 
highest redshift bin one) gives: 

log MB " (max) = m max = 0.34(±0.02)z + 8.99(±0.03); (11) 

Me 

while the maximum Eddington ratio (~ 0.45) is consistent 
with no evolution with cosmic time: 

(max) = 10 (emax) = 0.005(±0.006)z + 0.45(±0.01). (12) 

X>Edd 

Assuming that the shapes of Mbh and Eddington ratio 
distributions do not change with redshift, which is suggested 
by the fact that the value of er m and er e are independent 
of z (see Table [1} , Eq. [11] also describes the slope of the 
z-dependence of the mean quasar BH mass, and not only 



of the maximum mass of the quasar populations. Similarly, 
the mean (and not only the maximum) Eddington ratio is 
constant with redshift. 

4.2 Dependence of the results on the tblr — A£a 
calibration 

We tested whether a variation of the luminosity exponent 
of the virial calibration may affect the relative evolution in 
m max and e max derived in this paper. The L — tblr relation 
assumed in Eq. Q] (McLure & Dunlop 2004) is quite steep, 
although still consistent with the canonical j"blr oc AL a 6 
with b — 0.5 which is often assumed for idealised photoioni- 
sation. In order to quantify the effects that this has on the 
relative evolution in the maximum mass and Eddington ra- 
tio, we reproduced the analysis assuming the exponent on 
the luminosity term is b — 0.5 or b — 0.4, both of which are 
consistent to within 2a with the McLure & Dunlop (2004) 
calibration in which b — 0.62 ± 0.14. 

Figure [5] (upper panel) shows that the smaller is the 
luminosity exponent in the virial calibration, the flatter is 
the dependence on redshift of the BH masses. The results 
are only slightly affected by the choice of the L — tblr rela- 
tion, since the z-evolution determined assuming an exponent 
of 0.5 is consistent within lcr with the previous determina- 
tion, obtained assuming the McLure & Dunlop (2004) virial 
calibration. In Table [2] we give the best linear fit to the max- 
imum mass as a function of redshift for b = 0.5, b = 0.4 and, 
for comparison, b — 0.62. 

On the other hand, as regards the dependence on red- 
shift of the Eddington ratio, the picture is more delicate. 
This parameter appears to increase significantly with z as- 
suming a flatter L — tblr relation, while it was found to be 
constant with redshift within the assumed virial calibration 
(Eq. [T] see Fig. [5] lower panel) . In Table [5] the parameters 
of the best linear fit to the maximum Eddington ratio as a 
function of redshift are given for various values of the lu- 
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Table 2. Best linear fit to the maximum mass and to the maximum Eddington ratio as a function of redshift for various values of the 
luminosity exponent in Eq. [l](6 = 0.62, b = 0.5 and b = 0.4). For each value of b, the average value over all the redshift bins of the xi 
of the best fit probability density (Eq. [4} is also given. 
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Figure 8. Upper panel: small dots are the the virial BH masses 
of the Mgn sample given by Shen et al. (2008); the dash-dotted 
line reports the corresponding mean values. Red circles are our 
estimates of log — gH(gjfO an( j ^ ne re( j so \^ ij ne ; s ^he best fit 
reported in Eq. 1111 The Mbh vs. z dependence proposed by Fine 
et al. (2006) is the dashed line. Lower panel: small dots are the 
Eddington ratios for each source of the MgH sample; the dash- 
dotted line shows the corresponding average. Red circles are the 
maximum Eddington ratios and the red solid line is the best linear 
fit reported in Eq. 1121 



Figure 9. BH maximum masses (upper panel) and maximum 
Eddington ratios (lower panel) as a function of redshift for differ- 
ent values of the luminosity exponent in Eq. [T] 

best fit probability density are given for b — 0.62, b = 0.5 
and b = 0.4. It is apparent that the xt IS minimum for 
b — 0.62, giving a circumstantial independent support to 
the index proposed by McLure & Dunlop (2004, see Eq. [1} 
and, hence, to Eqs. Illl and 1121 



minosity exponent (b — 0.5, b = 0.4 and, for comparison, 
b = 0.62). 

Note that the assumption of a flatter L — tblr rela- 
tion leads to an increase of the residuals between the best 
fit probability density (Eq. [4| and the observed quasar dis- 
tribution. In Table [2] the xt values averaged over z of the 



4.3 Comparison with previous results 

We now compare our results with those obtained by McLure 
& Dunlop (2004), Fine et al. (2006) and Shen et al. (2008), 
focusing just on the slope of the Mbh and Eddington ratio 
evolution. 
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Fine et al. (2006), in order to reduce the effects of the 
Malmquist bias, concentrated on a subsample of the 2dF 
quasars with luminosity around L*(z) at each redshift. They 
observe a significant dependence of the quasar BH mass on 
redshift (M B h oc (1 + z) 3,3±11 ), but conclude that their 
result cannot directly be interpreted as evidence for anti- 
hierarchical "downsizing" because the ^-dependence they 
found is strongly dominated by the dependence on redshift 
of L* . For comparison, we repeated the entire fit procedure 
described above imposing that M BH ( m ax)(z) varied as pro- 
posed by Fine et al. (2006). Figure [7] shows the relative Xv 
values in each redshift bin: the fit appears inadequate if we 
assume their results. Note however that the error given for 
the evolution of the average BH mass of QSOs by Fine et 
al. (2004) is quite large, so that their results are consistent 
with those given here within la. 

McLure & Dunlop (2004) proposed that the observed 
increase of the quasar BH mass with redshift is entirely as 
expected due to the effective flux limit of the sample. To 
further test the possibility that the mean BH mass is in- 
dependent of redshift, we repeated again the fit procedure 
described above assuming that M BH (max) is constant over 
all the redshift bins. In Figure[7]we plot the relative xt, and 
again the fit is inconsistent with the data, giving further to 
evidence for an evolution of quasar populations with z. 

McLure & Dunlop (2004) and Shen et al. (2008), study- 
ing the SDSS DRI and DR5 samples, found that there is a 
clear upper mass limit of ~ 1O 1O M0 for active BHs at z > 2, 
decreasing at lower redshift. This trend is in good agreement 
with our results and can be explained assuming that the 
quasar number density peaks at a certain z pea k ~ 2 — 2.5 
and then flattens out (see for example Richards et al. 2006) . 
Around z p0 ak, both the high and the low mass end of the 
quasar BH mass distribution are more populated, so that the 
observation of very massive objects is likely (while low mass 
quasars cannot be observed due to the instrumental flux 
limit). Therefore, the slope of the "unbiassed" dependence 
on redshift of the maximum BH mass is raised below z pC ak 
and it is flattened above. This effect translates in evidence 
for a limiting BH mass for active BHs at z > 2, decreasing 
at lower redshift, that is apparent in all large samples of 
quasars. 

McLure & Dunlop (2004) observed a substantial in- 
crease of Eddington ratios with redshift and a similar trend 
is apparent from the sample of 2dF L* quasars of Fine et al. 
(2006) after correcting their data for the offset between the 
Mgn and Civ virial mass calibrations (see for example Shen 
et al. 2008). We suggest that the observed z-dependence of 
Eddington ratios is spurious, and that it is entirely domi- 
nated by the dependence on redshift of the average quasar 
luminosity due to the Malmquist bias. 

4.4 Discussion of the results 

Studying a sample of ~ 50000 SDSS quasars with 0.35 < 
z < 2.25, we obtained that the maximum mass of the quasar 
populations increases with z, while the maximum Eddington 
ratio is practically independent of redshift. 

These results are unaffected by the Malmquist bias and 
may be interpreted as evidence for evolution of the active BH 
population with redshift. Quasar samples at lower redshift 
are increasingly dominated by lower mass BHs, i.e. most 



massive BHs start quasar activity before less massive ones. 
This is indicative of anti-hierarchical "downsizing" of active 
BHs and it is in agreement with recent theoretical predic- 
tions by e.g. Merloni, Rudnik & Di Matteo (2006). 

Our findings may have implications for the study of the 
co-evolution of supermassive BHs and their host galaxies, 
even if they cannot be directly interpreted as evidence for 
evolution of the Mbh — Mbuige scale relation. There is obser- 
vational evidence that quasar host galaxies are already fully 
formed massive ellipticals at z ~ 2.5 and then passively 
fade in luminosity to the present epoch (e.g. Kotilainen et 
al. 2007, 2009; Falomo et al. 2008). Within this scenario, our 
results can be interpreted as an evolution with redshift of 
the parameter V = Msri/Mbuigei which would be 4-5 times 
larger at z ~ 2 than today. 

This is in good agreement with the results of Peng et 
al. (2006), who found that F is ~ 4 times larger at z ~ 
1.7 than today in a sample of 11 lensed quasar hosts. Our 
results are also consistent with Salviander et al. (2007), who 
examined a sample of SDSS quasars finding that galaxies of 
a given dispersion at z ~ 1 have BH masses that are larger 
by AlogAfaH ~ 0.2 than at z ~ (see Lauer et al. 2007 for 
a detailed discussion on the selection bias which may affect 
these results). 



5 SUMMARY AND CONCLUSIONS 

Starting from the SDSS DR5 quasar catalogue, we focused 
on the ~ 50000 objects for which Mgn line widths and 3000A 
monochromatic luminosities were available. This sample 
(0.35 < z < 2.25) was divided in 8 redshift bins. In each bin, 
the object distribution in the FWHM-luminosity plane was 
described in terms of a minimum luminosity limit (due to 
the instrumental flux limit), a maximum mass and a maxi- 
mum Eddington ratio. The assumed probability density was 
compared to the observed distribution of objects in order 
to determine the free parameters with a best fit procedure 
in each redshift bin. Errors on the best fit parameters were 
determined with Monte Carlo simulations. 

We tested the robustness of the procedure through some 
simulations, and showed that the maximum mass and the 
maximum Eddington ratio determined in each redshift bin 
depend neither on the quasar number density nor on the 
survey flux limit (which is responsible of giving raise to a 
Malmquist-type bias in the observed dependence on redshift 
of the mean quasar BH masses). 

We then studied the dependence on redshift of the max- 
imum quasar BH mass and of the maximum Eddington ratio 
and found clear evidence for evolution of the active BH pop- 
ulation with redshift. Over the redshift range studied, we 
obtained that the maximum mass of the quasar population 
depends on redshift as log(M BH (maxj/M©) = 0.34z + 8.99, 
while the maximum Eddington ratio is found to be practi- 
cally independent of redshift. 

This means that QSO samples at lower redshift are in- 
creasingly dominated by lower mass BHs, i.e. the more mas- 
sive a BH is, the earlier it starts quasar activity. Within 
a scenario in which quasar host galaxies are already fully 
formed massive ellipticals at z ~ 2.5, our results can be also 
interpreted as an evolution with redshift of the parameter 
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r = MBH/Mbuigc, which would be 4-5 times larger at z ~ 2 
than today. 
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